The Y(3940), Z(3930) and the X(4160) as dynamically 
generated resonances from the vector- vector interaction 



R. Molina and E. Oset 

Departamento de Fisica Tedrica and IFIC, Centra Mixta Universidad de Valencia- CSIC, Institutos de 
Investigation de Paterna, Aptdo. 22085, 46071 Valencia, Spain 

Abstract. We study the vector-vector interaction within the framework of the hidden gauge formalism for the sector with 
quantum numbers charm C = and strangeness S = in the region around 4000 MeV. We get five poles, three of which could 
be identified with the Y(3940), Z(3930) and X(4160). These poles appear with quantum numbers / = and J PC = 0++,2 ++ 
and 2 ++ , respectively, and can be considered as hadronic molecules made of D*D* , D*D*. 
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INTRODUCTION 

Recently, the B-factories at SLAC, KEK and CESR, which were originally constructed to test matter-antimatter 
asymmetries or CP violation, have discovered new hidden-charm states around the energy region of 4000 MeV. These 
new states do not seem to have a simple cc structure. They are naively called as XYZ particles. Some of them, which 
we consider in this manuscript, are the X(3940), the Y(3940), the X(4160) and the Z(3930). 

The X(3940) and the X(4160) were observed by the Belle Collaboration in the e + e~ -t J/y/X reaction as a DD* 
and D*D* mass peak respectively [1]. The X(3940) — > DD* has been observed but there is no signal for the DD or 
the ft)// y/ decays. Because the fact that the Tj r (15) and ri c (2S) were also produced in double-charm production, it was 
believed that the X(3940) could have J PC = 0~+, being a 3% charmonium state (tj"), but then the X(3940) should 
have a mass ~ 4050 MeV or even higher [2]. Concerning to the X(4160), the known charmonium states seen from 
a e + e~ — > J /y/D*D* reaction have J — 0, thus, if this state is identified with a 3 1 Sq (?]") or 4 1 Sq (?]"') charmonium 
state, the mass predicted in the first case will be smaller (4050 MeV) and higher (4400 MeV) in the second case [2]. 
On the other hand, the Y(3940) has been firstly observed by the Belle Collaboration and Babar has confirmed it in 
B —> KY — > KcoJ/y decays [3]. Nevertheless, the values for the mass and width reported by Babar are smaller than 
the Belle's values. Also, very recently, Belle has reported a narrow peak in the cross section for yy ^ ft)// y/ [4] that is 
consistent with the mass and width reported for the Y(3940) by the Babar group. This peak, dubbed as X(3915), has a 
mass M = 3914±4±2MeV and width r = 28± 12+| MeV. The measure of the mass favor the possible assignment 
of this new peak to the Y(3940) instead of the Z(3930). The production modes ensure that the X(4160) as well as the 
X(3940) and Y(3940) have C-parity positive. 

Recently, the CDF Collaboration at Fermilab has announced a narrow peak near the threshold, which is 

designated as Y(4140), observed in the B + -> J/\if<j)K + decay, it has a mass M = 4143 ± 2.9(stat) ± l.2(syst) MeV 
and a width T = 11 J^f^stat) ± 3.7 (syst ) MeV [5]. The width observed is quite different from the one reported by 
Belle for the X(4160) [6], suggesting that one is taking about a different state. 

Finally, the Z(3930) has been seen by Belle as a peak in the spectrum of DD mesons produced in yy collisions. The 
Belle measurements favors the 2 ++ hypothesis, making the assignment of the Z(3930) to the ^P^ix'cj) charmonium 
state possible [7, 8]. 

In [9, 10, 11], previous theoretical work about the origin of these particles has been done, supporting the interpreta- 
tion of the Y(3940) and the Y(4140) as D*D* and D*D* molecules respectively. 

In this paper we will propose a theoretical explanation on the nature of some of these XYZ states, providing structure 
and quantum numbers for them. Our work is based on the hidden gauge symmetry (HGS) formalism for the interaction 
of vector mesons, which was introduced by Bando-Kugo-Yamawaki [12]. Concretely, we study the case of a system 
of two vector meson with charm C — and strangeness 5 = around 4000 MeV (hidden charm sector). 



FORMALISM: THE VV INTERACTION 



Within the theorical framework, there are two main ingredients: first, we take the Lagrangians for the interaction 
of vector mesons among themselves, that come from the hidden gauge formalism of Bando-Kugo-Yamawaki [12]. 
Second, we introduce the potential V obtained from these Lagrangians (projected in s-wave, spin and isospin) in the 
Bethe Salpeter equation: 

T = (1-VG) _1 V , (1) 

where G is the loop function. Therefore, we are summing all the diagrams containing zero, one, two... loops implicit 
in the Bethe Salpeter equation. Finally, we look for poles of the unitary T matrix in the second Riemann sheet. All this 
procedure is well explained in [13, 14]. In these works, two main vertices are taken into account for the computation 
of the potential V: the four-vector-contact term and the three-vector vertex, which are provided respectively from the 
Lagrangians: 
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&m = y (VmVvV 1 V v - V V V^V^V v ) , (2) 

-4f ] = igifyVv - d v V^V v ) . (3) 

By means Eq. (3), the vector exchange diagrams in Fig. 1 are calculated. The diagram in Fig. Id) leads to a repulsive p- 
wave interaction for equal masses of the vectors [13] and only to a minor component of s-wave in the case of different 
masses [14]. Thus, the four-vector-contact term in Fig. la) and the t(u)-channel vector exchange diagrams in Fig. lc) 
are responsible for the generation of resonances or bound states if the interaction is strong enough. 
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FIGURE 1. Terms of the J£m Lagrangian: a) four vector contact term, Eq. (2); b) three-vector interaction, Eq. (3); c) / and u 
channels from vector exchange; d) s channel for vector exchange. 

As we are interested in the region close to the two vector meson threshold, the three momenta of the external 
particles can be neglected compared with the mass of the vector meson, and thus, we can make \q\/M ~ for external 
particles, which considerably simplify the calculation. 

In order to consider the pseudoscalar-pseudoscalar decay mode, and thus compare properly with the available 
experimental data for the width, a box diagram is also included, which we show in Fig. 2 for the particular case 
of the pp system. Crossed box diagrams and box diagrams involving anomalous couplings were also calculated in 
[13], but they were much smaller, specially in the case of the anomalous coupling, than the contributions coming from 
the box diagrams in Fig. 2, and they were not considered in later works [14, 15, 16]. In [13], two states were obtained 
within this formalism for the particular case of the pp system that were identified with the /2(1270) and the /o(1370) 
as pp molecular states. Whereas in [14], the procedure of [13] is extended to SU(3) and eleven states are found, five 
of them can be associated with data in the PDG: the /o(1370), /o(1710), / 2 (1270), $(1525) and ££(1430). In these 
works, the integral of two-vector-meson loop function can be calculated by means of the dimensional regularization 
method or by means of the use of a cutoff. This requires the introduction of one parameter, which is the subtraction 
constant ( or a cutoff parameter), that is tuned to reproduce the mass of the tensor states, whereas the other states are 
predicted. The box diagrams need two parameters, the cutoff for the integral, A, and one parameter, A^, involved in 
the form factor used. These parameters were considered to be around 1 GeV and 1.4 GeV respectively in [13] to get 
reasonable values of the /o(1370) and /2(1270) widths. In the later work of [14], these values also provide a good 
description of the widths for the other states. 

In [15], the work of [13, 14] is extended to SU(4) in the study of the p(co)D* system. Here, the V matrix is 
straigthforward extended to SU(4), as it was done in the work of [17]: 
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where the ideal mixing has been taken for ©, and 7/y/. The SU(4) symmetry is strongly broken by the explicit use 
of the masses of the vector mesons in the vector-exchange diagrams, where the terms of the potential that comes from 

the diagrams that exchange a D* meson are proportional to k = —f- ~ 0.15, and this gives rise to corrections of the 

m D* 

order of 10% of the p-exchange terms. In this work, three states are obtained with / = 1/2, 7 = 0, 1 and 2 respectively, 
whereas the interaction is repulsive for 1 = 3/2 and thus no exotics states appear. The subtraction constant a in the 
function loop is fine-tuned around its natural value of —2 [18, 19], in order to get the position of the ^(2460) state 
at the PDG. The interesting thing comes from the calculation of the box diagram, this box diagram only has 7 = 
and 2, therefore, the state found with 7=1 cannot decay to kD by means of this mechanism, which is the reason 
why we get a small width of 4 MeV, compared with that of the state for 7 = 2, around 40 MeV [15]. Thus, we find 
a reasonable explanation for why the D*(2640) has a small width, T < 15 MeV, when we associate to this state the 
quantum numbers J p = 1 + . The state for 7 = with mass around 2600 MeV is a prediction of the model and there is 
no data over this state up to now. 

As one can see, the model in SU(3) provides a good description of the properties of many physical states. Also, in 
SU(4) we have seen that interesting new states appear. Thus, in the next section, we study the case of C = and 5 = 
around the region of 4000 MeV. 
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FIGURE 2. Box diagrams for the case of the pp interaction. 



THE XYZ PARTICLES 

In view of the results obtained in the cases of the vector-vector interaction in SU(3) and the p{(0) D* system 
in [15] we make an attempt to study now the sector with C = and S = 0. In the case of / = we have 10 
channels: D*D*(4017), D*D*{4225), K*K*{1783), pp(1551), ©©(1565), 00(2039), J/yJ/y{6194), ©7/V(3880), 
07/i/a(4116) and ©0(1802),' and 6 channels for / = 1: D*/>(4017), K*K* {17 83), pp(1551), p©(1558), p7i/A(3872) 
and p 0(1795). The procedure is identical to the one explained in the previou section, see [16] for practical details, 
except by the complication added due to the many channels involved. The function loop is computed by means of 
dimensional regularization and we fix the value of p = 1000 MeV as in the work of [14]. We establish two different 
subtraction constants, one for the channels involving two light vector meson, which is called as 0£l, and the other 
one, for the channels involving two heavy vector meson: an- When one channel involves both kind of mesons, one 
heavy and one light vector meson, we also put o,l- The value of is set to get the position of the pole /2(1270) 
as OLl = —1.65, like it was done in the work of [14]. Thus, we get the correct position of the poles found in this 
work. The an constant is set to —2.07 in order to get the position of the pole found in S = around 3940 MeV. 
By the fact of having several channels, ones involving heavy mesons and the other light mesons, the g parameter in 
the Lagrangian must be different for the different cases. Thus, we put g = M p /{2f K ) = 4.17 for light mesons, and 
go = M D ,/{2f D ) = 6.9, g Ds = M D */{2f Ds ) = 5.47, gr]c = M J/¥ /{2f r]c ) = 5.2, when one D* , D* or 7/y/ meson is 
involved. The DD decay mode is also considered by means of a box diagram as it was done in [15], but the effect is 
small compared with the (light) vector-(light or heavy) vector decay modes. The results for the pole position obtained 
together with a possible assignment with some of the XYZ particles observed in the energy region of 4000 MeV is 
given in Table 1 . In order to do a proper comparison with the experiment, we give in Table 2 the XYZ discovered up 
to now in which we are interested, which is taken from [7]. 

As it is shown in Table 1, we find one pole with a mass ~ 3940 MeV for each spin 7 = 0, 1 and 2 and 1 = 0. The 
modules of the coupling contants to the different channels are obtained from the residues of the amplitudes and they 
are given in Tables 3 and 4. From Table 3, we see that these three poles around 3940 MeV are about 50 — 70% D*D* 
and 20 — 30% D*D*, whereas all the other channels account only for 10 — 30% or less. In Table 2 (experimental data) 
there are also three states around this mass. The state that we found for / = 1 is discarded to be associated with one 
of these states, because in our model this state decays to DD and not to ©7/y/, thus cannot be associated with the 
X(3940), since the DD decay has not been observed [1], whereas it cannot be associated to the Y(3940) either, because 



TABLE 1. Comparison of the mass, width and quantum numbers with the ex- 
periment. All the quantities are in units of MeV. 
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L J 


Theory 




Experiment 






0+[0++] 

o-[i+-] 

0+[2++] 
0+[2++] 


(Mass, Width) 
(3943, 17) 

(3945, 0) 
(3922, 55) 
(4157, 102) 


Name 
7(3940) 

%(3945)" 
Z(3930) 
X(4160) 


Mass 
3943 ±17 
3914.3+^ 
3914±4±2 : 

3929 ±5 
4156±29 


Width 

87 ±34 
33+^ 
28± 12+| 

29 ±10 
139111 3 


J PC 
J p + 

2++ 
J p + 


1"[2++] 


(3912, 120) 


%(3912)" 









TABLE 2. 


Properties of the candidate XYZ mesons. 






State 


M (MeV) 


T (MeV) 


J PC Decay modes 


Production modes 


Z(3930) 
X(3940) 
K(3940) 

X(4160) 


3929 ±5 
3942 ±9 

3943 ±17 
3914.3j^;J 

4156 ±29 


29 ±10 
37 ±17 
87 ±34 

33+f 
139+^ 


2++ 
J p + 
J p+ 

J p + 


DD 
DD* 
aJ/ij/ 

D*D* 


n 

e+e- ^y/v/X(3940) 
fl -> £7(3940) 

e+e- (4160) 



the experimental state decays to (oJ / y/. On the other hand, we find a width so large that it is likely to be assigned to the 
Z(3930). In that way, the state found for / = 1 is a prediction of the model. In Table 1 we call it as Y p (3912), where 'p' 
stands for prediction. For this state, we find that the D*D* channel contributes 60%, whereas the sum of the modules 
of the couplings for the rest of the channels gives <~ 40%. 

Although, at first one could think that the three states that we found around 3940 MeV for / = could correspond 
to those in Table 2, this assignment can lead to error. The reason is that all the three experimental states have C- 
parity positive, while our state for J = 1 has C-parity negative. From the experimental point of view, this is clear by 
the production mechanism: e + e~ has C-parity negative (it comes from a photon) and as the J/y has also C-parity 
negative, the X(3940) must have C-parity positive. Thus, the state that we found for / = 1 is also a prediction. 

In Table 2, we find that there are two different experimental measures for the mass and width of the Y(3940). One is 
done by Babar and the other one by Belle with bigger values for both magnitudes. The point in favor to the assignment 
of the state that we found with J PC = 0++ to the Y(3940) is the calculation of the r((3943,0+[0++]) ->• coJ/yr) that 
can be done straightforwardly by means of the formula: 

r((3943,0+[0 ++ ]) -+ coJ/w) = (5) 

where p is the momentum of co in the resonance rest frame. Taking the modules of the coupling \gYmj/ y I = 1445 MeV 
from Table 3, we obtain r((3943,0 + [0 ++ ]) -> coJ /y/) = 1.52 MeV, compatible with the expected experimental value 
for this decay, T(Y (3940) — > (Oj/y) > 1 MeV [7]. Then, we have found a natural explanation on why this rate is much 
larger than it would be if it corresponded to a hadronic transition between charmonium states. 

On the other hand, the two states that we find in / = for J = 2, with masses M = 3922 and 4157 MeV, are assigned 
to the Z(3930) and X(4160) respectively, by the proximity of the mass, width and quantum numbers. The structure 
of the state found with M = 4157 MeV is radically different to the other states. We see from Table 3, that the D*D* 
coupling accounts for the 49%, whereas the D*D* coupling gives only the 3%, being the sum of the contributions of 
all the other channels of <~ 48%. 



TABLE 3. Quantum numbers, pole positions and modules of the couplings |g;| in units of MeV for / = 0. 
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COifl 


0+[0++] 


3943 + ;7.4 
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8645 


15 


52 


1368 


1011 


422 


1445 


910 


240 
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3945 + i0 


18489 


8763 


40 
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3922 + j26 


21177 


6990 


44 


84 


2397 


1999 


1794 


2433 


3061 


789 


0+[2++] 


4169 + ;66 


1319 


19717 


87 


73 


2441 


3130 


2841 


2885 


5778 


1828 



TABLE 4. Quantum numbers, pole position and modules of the couplings 
| in units of MeV for / = 1. 
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PP 
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pJ/v P<P 


1"[2++] 


3919 + «74 


20869 


152 





3656 


6338 2731 



CONCLUSIONS 

We have studied the C — and 5 = sector around the energy region of 4000 MeV within the hidden gauge formalism 
combined with coupled channel unitarity. We find one state with mass ~ 3940 MeV for each spin / = 0, 1 and 2. Only 
the states with 7 = and 2 can be associated the XYZ particles with mass ~ 3940 MeV in the PDG. Concretely, we 
associate the J PC = 0++ state with the Y(3940), and the 2++ state with Z(3930). Nevertheless, the J PC = l+~ cannot 
be associated with known states and it becomes a prediction. The only state that we find for / = 1, with J PC = 2 ++ , 
cannot be associated either, because of the decay modes and the large width. This new prediction could be observed 
in the pJ/\j/ channel. For those states with masses around 3940 MeV, we find that the D*D* channel accounts for the 
50 — 70%. Our model predicts another state with mass around 4160 MeV and quantum numbers I G [J PC ] = + [2 ++ ], 
which we identify as the X(4160) by the proximity of the mass, width and C-parity. This state has a different structure, 
being the D * S D* channel the one contributing most to the generation of the state. 

The region around 3940 MeV is very interesting and there could be more resonances not yet seen in this region. 
The findings of this work should motivate the experimentalist to look into this region in the channels that involve light 
vector - light vector or light vector - heavy vector like K*K* and pj / y/. 
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